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Abstract 

Three-pion interferometry is investigated for new information on the space- 
time structure of the pion source created in ultra-relativistic heavy-ion col- 
Usions. The two- and three-pion correlations are numerically computed for 
incoherent source functions based on the Bjorken hydrodynamical model, over 
a wide range of the kinematic variables. New information provided by three- 
pion interferometry, different from that provided by two-pion interferometry, 
should appear in the phases of the Fourier transform of the source function. 
Variables are identified that would be sensitive to the phases and suitable 
for observation. For a positive, chaotic source function, however, a variation 
of the three-pion phase is found to be difficult to extract from experiments. 
Effects of asymmetry of the source function are also examined. 

PACS number(s): 25.75 Gz 
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I. INTRODUCTION 



Two-pion interferometry has been regarded as an important means of obtaining the 
space-time structure of dynamics involved in relativistic heavy-ion colhsions. Over the years, 
extensive studies of two-pion interferometry have been carried out theoretically and exper- 
imentally to investigate how much information the correlations of two emitted pions can 
provide. As experiments become refined, measurements of three-pion correlations should 
become feasible, hopefully providing new information. The first result of such measure- 
ments has been recently reported from the CERN NA44 experiment 0. 

In the last few years, some theoretical investigations have been made regarding three-pion 
interferometry 0,^-^. Though the treatments of three-pion interferometry and the issues 
involved are similar to those regarding two-pion interferometry, three-pion interferometry 
is technically far more complicated than the two-pion case due to the involvement of an 
additional momentum and also to new aspects of the particle correlations. Consequently, 
theoretical work so far has been on limited aspects of the interferometry, focusing on a few 
kinematic variables, such as the sum of three relative invariant momenta, over small ranges 
of their values. 

In the coming years, especially when RHIC becomes in operation, we expect that exper- 
iments on three-pion interferometry will become more detailed and will be made over larger 
ranges of various kinematic variables. We report here an investigation of three-pion interfer- 
ometry over a wide range of kinematic variables. Our major objective is to clarify how much 
new information we can extract from three-pion interferometry, regarding the space-time 
structure of the pion source created in ultra-relativistic heavy-ion collisions. This work is 
similar to that of Heinz and Zhang ||^ in its objective, but it differs in scope. We numerically 
calculate the two- and three-pion correlations for the same source functions and compare 
the two types of correlations. The calculations are carried out with various model source 
functions that are based on the Bjorken hydrodynamical model. Throughout this work, we 
assume pion emission from the source to be completely chaotic. We also neglect possible 
final-state interactions (including the Coulomb interactions) between the emitted pions and 
the source. 

The new aspect of the correlations that three-pion interferometry can provide is the phase 
of the source function's Fourier transform. The information content of the phase differs from 
that of the magnitude of the Fourier transform. Since two-pion interferometry can provide 
only the magnitude, we hope that three-pion interferometry can provide new information. 
To this end, we carefully identify the variables that are sensitive to the phase and suitable 
for observation. Since the phase is expected to be greatly affected by asymmetry of the 
source function, we also examine that issue. 

In Sec. II, we summarize the formulation of both the two- and the three-pion inter- 
ferometries, mainly to define the correlation functions used. In Sec. Ill, we describe the 
various source functions that are used in the calculation, based on the Bjorken hydrodynam- 
ical model. The choice of optimal variables in the three-pion correlations is given in Sec. 
IV, and various numerical results are described in Sec. V. Discussion and conclusions are 
presented in Sec. VI. 
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II. FORMALISM 



The two- and three-particle correlation functions that we discuss are reasonably well- 
known, but for clarity we sketch the formalism of, in our case, bosons, and we define the 
correlation functions and various variables that appear in this work. 

The field operator ip{x) of emitted bosons, such as pions and kaons, obeys the Klein- 
Gordon equation, 

(9' + m2)^(x) = J(x), (1) 

where m is the mass of the particle, and J{x) the source current |Q. The incoming and 
outgoing states of ip{x) are specified in terms of the creation and annihilation operators 
and a for these states. They are related to J{x) as 

aout(p) = ain(p) + i f (fx ^ =J(x)e~^P-^, (2) 
J v'(27r)3 . 2up 

where po is on-shell and uOp = v^p^ + m^. 

We define an important quantity in this work, a source function, S{x,K), in terms of 
J(x): 

Six, K)= j ( Jt(x + |) J(a; - |))e^^-^A, (3) 

where (...) denotes the quantum ensemble average over various incoming states. The Fourier 
transform of S{x, K) is defined as 

p(g, K) = j S{x, K)e''^-''d^x (4) 
^fiq,K)exp{t<Piq,K)). (5) 

where f{q, K) is the amplitude and K) the phase. For a real S{x, K), we have /(g, K) = 
f{—q,K) and (j){q,K) = —(f){—q,K), and thus (f){0,K) = 0. Boson spectra and correlation 
functions will be expressed in terms of p{q,K) or S{x,K), and (f){q,K) is written in the 
small-momentum expansion as 

K) = {x)k -q-liiq- ^k?) ^ + O(g^), (6) 

where {...)k = / S{x, K)...d'^x, and xk = x — {x)k- Since (f){q,K) is written in terms of 
odd moments of x, an asymmetric (about x) source causes a strong, nonlinear dependence 
of (j){q, K) on q. In this work, we examine the S{x, K) of simple forms as discussed in Sec. 
III. 

The one-particle spectrum is given by 

Wi{pi) = (aLt(Pi)aout(pi)) 

= (2^0W/^^"'^^^^'" 
1 



Higher-order correlations of J(x) are assumed to satisfy Gaussian reduction, such as 

(jt(xi)jt(x2)J(yi)J(l/2)) = 

The two-particle spectrum is then written as 



(8) 



W^2(pi,P2) = (aLt(Pi)«out(pi)aLt(P2)aout(P2)) 

= Vri(pi)W^l(p2) ^ 

= W^i(pi)Vri(p2) 



(27r)6 . 2cup, ■ 2cup, 
1 



(9) 



(27r)6 . 2cup, ■ 2ujp^ 

where for the pair of the i-th and j'-th momenta, we define their average and relative momenta 



as 



2 



Qij Pi Pj y 

respectively. The two-particle correlation function is then expressed as 

W^2(Pl,P2) 



(10) 



C2(Pl,P2) 



W^l(Pl)W^l(P2) 
/l2 



1 + 



fuf'. 



22 



(11) 



Hereafter, we denote fij for /(%■, Kij) and for (f){qij, Kij). Equation (|TI|) shows that the 
two-particle correlation function is independent of (pij. 

Similarly, the three-particle spectrum and correlation are given by 



W^3(Pi, P2, Ps) = (4ut(Pi)aout(Pi)aLt(P2)aout(P2)aLt(P3)aout(P3)), 



(12) 



and 



C3(Pl,P2,P3 



W^3(Pl,P2,P3 



W^l(Pl)W^l(P2)W^l(P3) 



1 



fl2 , /23 , f? 



+ 



where 



/11/22 722/33 
-F2F3Cos<l>, 



012 + 023 + 031 
/12/23/3I 



+ 



31 



/33/1I 



(13) 



/11/22/; 



33 



(14) 



From the preceding discussion of (pij, we have 
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^{Pl =P2=P3) =0 



(15) 



for a real source function. When {x)ki2 ~ {^)k23 ~ (3^)^-31, Eq. (6) yields 

$ ^ -^[((gi2 ■ S:kJ')^^^ + ((923 ■ S^kJ')^^^ + {{qsi ■ S^kJ')^J ; (16) 

$ is thus expected to vary prominently when the source is asymmetric (in x). 

Equation (|13D shows that the three-particle correlation function depends on cos$, which 
is absent in the two-particle correlation. The rest of the three-particle correlation is expressed 
in terms of /j^'s and can thus be determined from the two-particle correlations. 

The three-particle correlation is a function of three momenta, Pi,P2, and ps. For con- 
venience, we introduce the average total momentum, 

K = Pl±Pl±Pl. (17) 

Since qi/s satisfy the identity, 

qi2 + + qsi = 0, (18) 

we have three independent momentum variables, which we will take to be K, qi2, and ^23- 

Among various choices of three variables, another convenient choice is a set of -ft'j/s. For 
completeness, we show the relations between the first set and -R'j/s: 

K12 = K + ^{q23- qsi), K23 = K +^{q3i - qi2), ^31 = ^+^(^12-523), (19) 



satisfying 

K = (20) 



K12 + i^23 + ^31 



III. SOURCE FUNCTION 

We apply the Bjorken hydrodynamical model to describe the evolution of the hot region 
created in heavy-ion collisions. Based on the Cooper-Frye spectrum, one writes the source 
function as @] 

K-U\ 



S{x, K)d^x = F{t, ri)dTexp I —777^ \ x K ■ daexp 



F{T,ri)dT exp I x rriT cosh.{Y — rj)e ^"^^("^ ^\drjd^x.T, (21) 

\ Zltrp 



where r, rj, and xj- denote the coordinate variables defined as r = \Jt'^ — z^, V = \ log (|rf 
and = {x, y) , respectively; mrp and Y are the momentum variables defined as m-p = 
Kq — Kl and Y = \ log , respectively. Note that we specify the momentum 
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variables for a momentum pair, pi and pj, by the subscript ij, such as rriTij, and those for 
three momenta with no subscript such as rriT. We set the 2;-axis to be the beam direction, 
and the x-axis parallel to Kj^, or perpendicular to the beam axis. If^ is a four-velocity of flow, 
and is f/^ = (^,0,0, ^) in the Bjorken hydrodynamical model. The profile function F(r, 77) 
determines the source shape in the T-rj space. The shape along the transverse direction is 
taken to be of a Gaussian form. Tf is the freeze-out temperature, and da^ is the measure 
of freeze-out hypersurface. We assume that freeze-out occurs on the hypersurface where r 
is constant. 

In this work, we examine the following five forms of the profile function. 
Simple (box-type) profile: 



Gaussian profile: 



^^'^'''^"[0 (otherwise) ^^^^ 



F(r, r/) = 5(r - tq) exp (--^] (23) 



Heinz profile: 




(27r)Ar 



Exponential profile: 



Theta profile: 




(24) 



(26) 



IV. OPTIMAL VARIABLES 

The three-particle correlation function depends on three momenta, which have nine com- 
ponents. As noted in Sec. II, we choose the three momenta to be K, qi2, and 523- In this 
work, we focus on the dependence of the correlation functions on relative momenta of the 
emitted bosons by fixing the value of K. This leaves the two relative momenta to be the 
remaining variables. 

In order to identify new information in the correlation function, we should choose the 
variables that could provide the most rapid variation of cos Figure illustrates the 
variation of cos$ as a function of gfg Q23 "with all other components of the relative 
momenta set to be zero. Here, we use the Heinz profile with the parameters of tq = 6 fm. 
At = 1 fm, Af] = 1.2, and ttit = Tj = 140 MeV. Figure |I| shows that cos$ is unity along 
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the gf2" g|3-axes and also along the line of gfg = "Q'Isj the gf^^-axis (because of Eq. 
(16)). Figure |l| also shows that cos$ varies most prominently along the lines of gfg — ^23 ^ 
Q23 = '?3i5 q^i = qf2- Because of symmetry, however, it is sufficient to examine the 
variation around one of the lines, such as the line of — ^23- 

The reason why the line of — Q23 provides a prominent variation of cos $ is seen as 
follows. When qij <^ Kij and K12 ~ K23 ~ -ft'si, we have 

$ ~ 0(gi2) + 0(^23) + 0(-gi2 - 923), (27) 

where the (weak) dependence on Kij in cf) is not explicitly shown. Setting = Q'z cos 9 and 
?23 = '?2 sin 6*, we find 

— ^ - q^ sin 6'0'(gi2, ?i2, 5^ cos 6*) + cos Ocjy'iq^s, q^^ qz sin 6) 

- q,{cose - sme)(j)'{-q^2 - ^23' -?i2 " ^23, -^^(cose + sin^)), (28) 

where 0'(q) = dcp/dq^. When qf2 = ^23 112 = ?23) d^/dO vanishes at 6* = 7r/4, or 
I12 ~ ?23' yielding the minimum. It follows then that the optimum choice is qi2 = q23- 

With this constraint, the variables are now reduced to a single relative momentum. 

The source functions under consideration (as listed in Sec. Ill) are almost symmetric 
about the beam and transverse axes, but they can be asymmetric about the time-axis. 
Which component of gi2(= 523) should we choose so as to describe most effectively the 
asymmetric property of the source functions? g°2 is an obvious choice, but we find that gf2 
or qf2 is the most convenient variable. From Eq. (10), we have K12 ■ qu = 0, which gives 
gi2 ■ x = qi2iiKl2/ K^2 ~ ^) + ?i2(^-^i2/-^i2 — x) — quV- Thus cos$ as a function of gf2 
qi2 refiects time-axis asymmetry. Furthermore, for finite Y'12, our source functions are not 
simply a function of without symmetry along the longitudinal direction. Because of this, 
q^2 may be the variable more suitable for identifying the time-axis asymmetry of the source 
function. In the following, we will examine both g^2 q^2- 

qi2 is parallel to Kfg; while qi2 is perpendicular to it. Following the practice in the litera- 
ture on two-pion interferometry, we will also denote g^2 ^^"^ 112 ^ lout and gside, respectively. 
Note that when the dependence of cos $ on g^2 or is examined, the other components of 
qi2 are fixed in our calculation. 



V. NUMERICAL RESULTS 

Figures ^ ^, ^ |^, and [1^ illustrate the dependence of /12 and 0i2 on g^2 the Heinz, 
Simple, Gaussian, Theta, and Exponential profile functions, respectively. Each figure is 
shown for Y12 = 0, 1, and 2. In the figures, we set g°2* = ^12''^ = O5 ^ti2 = 140 MeV, Tf = 
140 MeV, Arj = 1.2, and Tq = 6 fm for all profile functions, and also At = 1 fm for the 
Heinz, Theta, and Exponential profile functions. These parameter values reasonably satisfy 
the recent CERN-SPS experiment on Pb + Pb at 158 GeV/A [§. Note that \fu\^ + 1 
describes the two-particle correlation since fu is normalized to be unity at gfg = 0) and that 
is not observed in the experiment. 

Figures H, |^, |^, ||, and |Tl| show cos $ and (half of) its coefficient in the three-particle 



correlation function, F3, as functions oi q^ { = qf2 = 123)- Each figure is shown for F = 0, 1, 
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and 2 and rriT = 140 MeV with the parameter values of Tf, At], tq, and At as in Figs. ^, 
I, I, and 0. 

In Figs. ^, 0, 1^, and |11|, we see that for all profile functions, cos $ = 1 at y = 
0, and cos $ becomes more prominent as Y increases, cos $ tends to be smaller for the 
asymmetric profile functions (such as the Simple and Exponential profile functions) than 
for the symmetric ones, though the difference between them is not substantial. $ tends to 
deviate from zero more slowly than the phases of Pi/s, (pij^s. For example, we see in Figs. 
4 and 5 that when (j)i2 reaches n/2 around gfj — 1^0 MeV, we still have cos$ > 0.5. The 
reason for this is general because $ is defined as $ = 012 + 023 + 031 and is constrained by the 
identity, gi2 + ^23 + Q'31 = 0. In fact, the small- momentum expansion of Eq. (|^) applied to $ 
yields that the term linear in q vanishes, as discussed in Sec. II. Furthermore, when $ starts 
to deviate from zero, its coefficient in the three-particle correlation function, 2F3, tends to 
become small. In fact, we see in the figures that 2F3 gets halved for all profile functions 
before cos $ decreases to 0.9, and even nearly vanishes when cos $ decreases further. We 
expect that in actual experiments it will be difficult to identify cos $ with a value much 
smaller than unity. This is the major finding of this work. 

We find the same difficulty when we choose gout as the independent variable. Figure ^ 
shows /12 and 0i2 as functions of g°2* Heinz profile for rnxu = 140, 200, and 300 

MeV, with tq = 6 fm, Y12 = 1.5, Arj = 1.2, At = 1 fm, and Rt = 5 fm. Figure 13 also shows 
cos$ and F3 as functions of gout(= Q'12* = ^23*) three values of rriT with Y = 1.5 and 

the same values of Tf, Ar], At, tq, and Rt- We see in Fig. |13| that when cos$ starts to 
deviate from unity, its coefficient, 2F3, becomes small quite rapidly, as in the previous cases 
of being the independent variable. Note that F3 decreases more quickly as a function of 
gout than of g^, as seen in Figs, [l^ and ^ The quick reduction occurs because of the factor 

exp (^—-^r^, when the transverse momentum, gout, is used as the variable. Rt = 5 fm is 

chosen to reproduce experimental results, but we find that even if a much smaller value of 
Rt is used, F3 becomes quite small when cos $ is off unity. 

The difficulty also remains at different values of tq. Figures n and ^ show /12, 0i2, 
cos$, and F3 as functions of gfa for (tq, Ar) = (6.5, 0.65), (4,8, 2.8), and (3.2, 3.2) fm, with 
rriT = mTi2 = 140 MeV, Tf = 140 MeV, Y = Y12 = 1.5, and Ar] = 1.2. These parameter 
sets yield the recently measured value of the longitudinal size, -Riong = 3.5 fm |^. Figure ^ 
shows that cos $ depends on tq rather strongly. But the variation of cos $ takes place where 
F3 is small, and will be difficult to observe. 



VI. DISCUSSIONS AND SUMMARY 

We have investigated the three-particle correlation function for chaotic source, using 
various profile functions of the source in the Bjorken hydrodynamical model. In all profile 
functions we have examined, the coefficient of cos $ in the three-particle correlation function, 
2F3, decreases quite rapidly, as cos$ decreases from unity. Extraction of cos$ is thus 
difficult when cos $ decreases from unity. This result is in agreement with Heinz and Zhang 
0, who previously made a similar investigation for small relative momenta. 

In order to clarify why cos $ is difficult to observe, we first identify the portion of the 
source function S{x) that most strongly influences the phase, 0. appears in the Fourier 
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transform of S{x) as p{q) = pR^q) + ipi{q) = fe^'^. As noted previously, is an odd function 
of q, and 0(g = 0) = for a real source function. The imaginary part of p{q), Pi{q), thus 
represents the behavior of relevant to the present discussion and is an odd function of q. 
The portion of S{x) that most strongly affects is thus its odd-function part, Sodd{x), in 

S{x) = Seven{x) + Sodd{x) with Sodd{-x) = -Sodd{x) and Seven{-X) = Seven{x). 

Consider that Seven{x) and Sodd{x) vary over distances of typical scale, 6xeven and Sxodd, 
respectively. Here, dxeven and 6xodd correspond to the spread widths of Sevenix) and the 
distance to maximum of Sodd{x), respectively, as Fig. |16| illustrates. pR^q) and pi{q) also 
vary over the scales of Sx~^^^ and Sx~]^^, respectively. As increases, pi{q) increases. The 
region where pi{q) is appreciable is g > Sx~^^. The amplitude / becomes small in the region 
where q is larger than 

Thus, if Sx~j^ < Sx~^^^, the variation of is not difficult to observe. This is, however, 
untenable as long as S{x) is positive everywhere. Consider Sodd{x) < in some region: we 
then have Seven{x) > \Sodd{x)\ in that region, and furthermore, Sxeven > ^Xodd- The last 
inequahty contradicts the desired inequality. 

In this discussion, we have assumed the source function, S{x,K), to be positive ev- 
erywhere, as done in common practice. How realistic is this assumption? If S{x, K) is a 
statistical phase-space distribution, S{x, K) should not be negative everywhere. Generally 
speaking, however, S{x, K) can be locally negative, provided its integrals over a; or i^T be 
positive. This property is similar to, for example, that of the Wigner function. Furthermore, 
S'(x, K) that is locally negative suggests that it would be associated with some dynamics 
involving quantum correlation. This aspect of the interferometry is currently under investi- 
gation. 

Measurement of the three-pion correlations from the CERN NA44 experiment has been 
recently reported |I| for the total relative momentum of up to about 300 MeV, and the 
extraction of cos $ up to about 60 MeV. This total relative momentum is too small to 
observe the variation of cos$ that is examined in this work. Taken as a constant, cos$ is 
found to be about 0.2 [|l| and is much smaller than the value of unity for a chaotic source, 
as discussed in this work. The reason for the small value is unknown, but it may be because 
of partial coherence of the source, or even because of a more exotic reason . 



In summary, we have investigated three-particle interferometry from chaotic sources us- 
ing the Bjorken hydrodynamical model. The optimal variables are identified suitable for 
extracting new information through the phase of the source function's Fourier transform. 
The three-particle correlations are calculated over a wide range of the kinematic variables. 
A variation of the three-pion phase is found to be difficult to observe experimentally be- 
cause its coefficient, as it appears in the three-particle correlation function, becomes small 
in that region. This is the case if the source function is positive everywhere as convention- 
ally assumed, and suggests the interesting possibility of the source function being locally 
negative. 
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FIGURES 




FIG. 1. cos$ as a function of Q23 fo^" Heinz profile function. The contour curves on 

the qf2 — Q23 plane are shown for the value of cos $ from to 1 by the step of 0.1. See the text for 
the parameter values used. 
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FIG. 2. /12 and (pi2 as functions of for Y12 = 0, 1, and 2 in the Heinz profile function. See 
the text for the parameter values used. 
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FIG. 3. cos$ and F3 as functions of Qz with y = 0, 1, and 2 for the Heinz profile function. See 
the text for the parameter values used. 




FIG. 4. /12 and ^12 as functions of grf2 with Y12 = 0, 1, and 2 for the Simple profile function. 
See the text for the parameter values used. 
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FIG. 5. cos$ and F3 as functions of Qz with y = 0, 1, and 2 in the Simple profile function. See 
the text for the parameter values used. 
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FIG. 6. /12 and 4>i2 as functions of qf2 with Y12 = 0, 1, and 2 in the Gaussian profile function. 
See the text for the parameter values used. 
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FIG. 7. cos$ and F3 as functions of Qz for y = 0, 1, and 2 in the Gaussian profile function. 
See the text for the parameter values used. 
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FIG. 8. /12 and 4)12 as functions of ql^ for 1^12 = 0, 1, and 2 in the Theta profile function. See 
the text for the parameter values used. 
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FIG. 9. cos$ and as functions of Qz for Y = 0,1, and 2 in the Theta profile function. See 
the text for the parameter values used. 
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FIG. 10. /i2 and 0i2 as functions of qf2 for Y12 = 0, 1, and 2 in the Exponential profile function. 
See the text for the parameter values used. 
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FIG. 11. cos $ and F3 as functions of Qz for y = 0, 1, and 2 in the Exponential profile function. 
See the text for the parameter values used. 
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FIG. 12. /12 and ^12 as functions of q^^'^ for rriTu = 140, 200, and 300 MeV in the Heinz profile 
function. See the text for the parameter values used. 
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FIG. 13. cos $ and F3 as functions of qout for uit = 140, 200, and 300 MeV in the Heinz profile 
function. See the text for the parameter values used. 
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FIG. 14. /12 and 0i2 as functions of for (-^o, Ar) = (6.5, 0.65), (4.8, 2.4), and (3.2,3.2) fm for 
the Heinz profile function, shown by solid, dashed, and dotted curves, respectively. See the text 
for the parameter values used. 
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FIG. 15. cos$ and F3 as functions oi for (rcAr) = (6.5, 0.65), (4.8, 2.4), and (3.2,3.2) fm 
for the Heinz profile function, shown by sohd, dashed, and dotted curves, respectively. See the text 
for the parameter values used. 




FIG. 16. Distances of typical scale over which Seven and Sodd vary. 
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